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FOREWORD 

This  technical  report  summarizes  research  performed  in-house  at  the 
High  Speed  Aero- Performance  Branch,  Aeromechanics  Division,  Flight 
Dynamics  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Wright- 
Patterson  Air  Force  Base,  Ohio.  The  work  was  performed  under  Project  No. 
2404  "Aeromechanics,"  Task  No.  240407  "Aeroperformance  and  Aeroheating 
Technology,"  Work  Unit  24040732  "Reentry  Vehicle  Performance  Analysis 
Methods."  The  study  period  was  March  1979  to  September  1981. 

The  report  was  written  by  David  T.  Johnson,  Charles  H.  Bursey,  and 
Richard  C.  Nash  of  the  Flight  Performance  Group.  Data  runs  were 
accomplished  with  the  assistance  of  Everett  M.  Marshall.  Special  thanks 
are  extended  to  Gerald  L.  Burke  for  the  heating  methods  of  Section  III 
and  to  Mr.  Fay  Young  of  ASD/AD  for  the  programming. 
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SECTION  I 
INTRODUCTION 

The  ability  of  high  speed  maneuvering  vehicles  to  perform  military 
missions  ncs  been  a  subject  of  interest  over  the  last  two  decades.  Studies 
have  been  made  of  concepts  ranging  from  space  transportation  and  on-demand 
rescue  to  boost  glide  maneuvering  missiles.  A  great  deal  of  effort  has 
gone  into  the  development  of  aerodynamic  shapes  for  these  purposes.  To 
determine  which  concepts  were  most  suited  some  sort  of  performance  measure 
was  needed.  With  the  advent  of  the  digital  computer  as  a  practical 
engineering  tool  a  performance  "footprint'1  as  a  measure  of  maneuvering 
capability  was  selected.  This  cross  range-down  range  footprint  of  a 
vehicle  concept  shows  the  landing  or  impact  versatility  which  can  be  gained 
by  maneuvering  in  the  earth's  atmosphere.  The  footprint  is  directly 
dependent  on  lift-to-drag  (L/D)  capability  of  the  vehicle  and  shrinks  to 
the  ballistic  impact  point  in  the  limit  where  L/D  =  0.  Generally  the 
first  of  these  footprints  to  be  studied  used  a  constant  L/D  and  bank  angle 
with  the  footprint  generated  by  varying  the  bank  angle  in  successive 
cases.  This  was  followed  by  a  staged  bank  angle,  i.e.  bank  set  at  a 
constant  value  until  the  heading  changed  90°  and  then  set  to  zero  during 
the  remaining  glide  range,  again  with  the  footprint  generated  by  varying 
the  constant  bank  portion  parametrically.  This  staged  bank  approach 
yielded  slightly  larger  cross-ranges  than  the  constant  bank.  Variable 
L/D  aerodynamics  were  also  generally  factored  into  the  analysis  at  this 
time  (CL  and  CQ  varying  as  a  function  of  altitude  and  Mach  number). 

The  next  step  in  footprint  development  evolved  with  the  construction 
of  flight  path  optimization  methods.  Computer  programs  were  generated 
and  applied  to  the  problem  of  maximizing  cross-range  (unconstrained). 

The  optimization  program  used  within  the  Flight  Dynamics  Laboratory  was 
based  on  the  steepest  descent  technique  and  is  still  available  although 
not  much  used  today.  Reference  1.  The  optimum  cross-range  solutions 
obtained  from  this  program  were  higher  than  any  developed  by  the  previous 
methods.  However,  to  obtain  a  solution  for  just  the  maximum  cross-range 
point  on  the  footprint  through  steepest  descent  optimization  was  time 
consuming  in  terms  of  computer  resources,  and  to  optimize  the  complete 
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footprint  for  each  vehicle  and  entry  condition  was  out  of  the  question. 
Footprints  were  hence  estimated  by  fairing  a  curve  through  the  maximum 
cross-range  point  and  the  parametric  points  generated  by  the  staged  bank 
procedure  described  previously.  Reference  2  contains  a  partial  list  of 
references  in  the  performance  area. 

The  objective  of  this  effort  was  to  develop  a  reliable  method  for 
rapidly  analyzing  entry  vehicle  performance  (see  also  Reference  3). 
Development  of  the  control  laws  stressed  the  use  of  realistic  vehicles 
in  the  medium  to  high  L/D  range,  that  is  those  which  are  supported  by 
analysis  and  wind  tunnel  data,  rather  than  hypothetical  vehicles  which 
were  much  used  in  the  past.  The  impact  of  heating  considerations  was 
also  addressed. 

The  approach  used  was  to  develop  optimal  cross-range  solutions  for  a 
footprint  both  without  constraints  and  with  down-range  constraints. 
Parameter  optimization  on  bank  angle  was  used  to  develop  these  solutions 
using  the  optimization  procedures  described  in  References  4  and  5  and 
the  trajectory  program  described  in  Reference  6.  Observation  of  the 
variation  of  the  optimal  bank  angle  control  schedule  with  respect  to 
turn  angle  then  led  to  the  development  of  a  simplified  control  for  rapid 
footprint  generation.  Comparisons  between  the  simplified  control  results 
and  the  optimization  results  were  then  made  for  three  conceptual  vehicles. 
Section  II  contains  these  results. 

In  Section  III  optimal  solutions  with  heating  constraints  were  sought 
and  the  effect  on  the  footprint  was  noted.  The  heating  method  employed 
was  a  reference  heat  transfer  procedure  similar  to  that  used  in 
Reference  7  in  conjunction  with  the  optimization  and  trajectory  programs. 
In  Section  IV  a  display  concept  is  presented  to  aid  in  translating 
footprint  data  to  functional  performance.  This  procedure  allows  plotting 
and  manipulating  of  data  on  world  maps  and  can  be  used  separately  or  in 
conjunction  with  the  performance  calculation  procedure.  Study  conclusions 
are  shown  in  Section  V.  Mathematical  details  are  shown  in  Appendix  A, 
vehicle  aerodynamics  in  Appendix  B,  and  a  footprint  data  base  in 
Appendix  C. 


2 


AFWAL-TR-82-3036 


SECTION  II 

FOOTPRINT  GENERATOR  DEVELOPMENT 

1.  REENTRY  VEHICLES 

Three  lifting  entry  vehicles  with  maximum  lift/drag  (L/D)  ratios  of 
1.89,  2.2,  and  2.6,  respectively,  were  selected  for  analysis.  These  L/D 
values  were  measured  at  Mach  =  20  and  200,000  feet  altitude.  All  three 
vehicles  are  realistic  in  terms  of  system  type  analysis  supported  by  a 
broad  base  of  aerodynamic  experimental  data.  Sketches  of  these  config¬ 
urations  are  shown  in  Figures  la  through  1c.  Aerodynamic  characteristics 
used  in  this  analysis  are  shown  in  Appendix  B.  These  data  were  obtained 
from  References  8  through  10.  The  Reference  9  data  were  supplemented  by 
theoretical  calculations  at  Mach  =  10  and  20.  A  wing  loading  of  70  psf 
was  used  for  the  low  L/D  vehicle  and  40  psf  was  used  for  the  other  two. 

2.  ASSUMPTIONS 

Atmospheric  entry  was  initiated  at  60  nautical  miles  (364,566  feet) 
altitude  heading  east  on  the  equator  (heading  a  =  90  degrees).  An  earth 
relative  velocity  of  24,480  feet/second  was  used.  This  corresponds  to  an 
inertial  velocity  of  26,032  feet/second,  which  is  slightly  above  circular 
velocity.  Although  aerodynamic  heating  considerations  limits  the  entry 
path  angle  corridor  that  is  practical  for  lifting  entry  vehicles,  entry 
paths  angles  from  -2.0  to  -5.0  degrees  were  investigated  for  the  purpose 
of  methodology  development.  An  entry  path  angle  of  -2.0  degrees  was 
used  for  the  purpose  of  footprint  generator  development  and  for  trajectori 
with  heating  constraints  imposed.  A  rotating  oblate  earth  with  the  1959 
standard  atmosphere  was  used  in  the  trajectory  computations.  Coordinated 
turns  (zero  side-slip  angle)  were  assumed.  Unconstrained  trajectories 
were  flown  at  zero  bank  angle  to  the  initial  pull-out  altitude  (y  =  0); 
bank  angle  was  then  allowed  to  vary  for  the  remainder  of  the  trajectory. 
Unconstrained  optimal  solutions  and  maximum  cross-range  under  down-range 
constraints  used  the  angl e-of-attack  for  maximum  L/D.  This  varied  with 
Mach  and  altitude  for  the  particular  vehicle.  The  optimal  solution  with 
heating  constraints  in  Section  III  varied  the  bank  angle  after  the  initial 
pull-out  and  the  angle-of-attack  over  the  entire  trajectory. 
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3.  EQUATIONS  OF  MOTION 

The  equations  of  motion  are  solved  in  an  earth-referenced  cartesian 
coordinate  system.  The  relation  between  local -geocentric-horizon, 
inertial  and  earth-referenced  coordinate  systems  is  illustrated  in  the 
following  sketch. 

Local -Geocentric 
Horizon  Coordinates 


Inertial  Coordinates 
Earth-referenced  coordinates 

Local -geocentric-hori zon  coordi nates 

The  origin  of  these  systems  lies  on  the  polar  axis  in  the  equatorial 

plane.  The  Z  and  axes  are  coll  inear  with  the  polar  axes  and  positive 

toward  the  south  pole.  The  Xg  axis  is  in  the  equatorial  plane  and  is 

fixed  at  the  longitude  of  the  vehicle  at  time  equal  zero.  The  Yg  axis  is 

positive  to  form  a  right-handed  system.  The  inertial  coordinates  X,Y,Z 

and  the  earth-referenced  coordinates  X„,  Y  .  Z  coincide  at  time  zero. 

e  e  e 

The  earth-referenced  axes  rotate  with  the  planet  whereas  the  inertial 
axes  remain  fixed. 


X,  Y,  Z, 
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The  equations  of  motion  solved  in  an  earth-referenced  cartesian 
coordinate  system  are  as  follows: 


X  =  F  /m  -  2  y  w  +  x  to  „ 
e  x  -'e  p  e  p 

e 


Y  =  F  /m  +  2  x  <jj  tyi) 
e  yg'  e  p  Je  p 


(1) 

(2) 


Ze  =  Fz'- 
e 


(3) 


These  equations  were  integrated  to  obtain  the  velocity  X  ,  Y  .  Z  and 

6  6  6 

position  X  ,  Y  ,  Z  . 

e  e  e 


The  transformations  necessary  to  relate  the  axes  systems  employed  are 
developed  in  Appendix  A. 


4.  OPTIMIZATION  PROCEDURE 

The  control  variables  used  for  optimization  were  bank  angle  (<J>B)  and 
angle-of-attack  (a).  These  variables  were  input  to  the  computer  program 
in  tabular  form  as  "control  parameters".  The  parameter  optimization 
procedures  then  vary  these  control  parameters  to  maximize  or  minimize  a 
given  function.  The  control  used  by  the  program  is  formed  by  a  linear 
interpolation  between  the  discrete  control  parameter  points  varied  by  the 
optimizer.  In  this  section,  eight  parameters  were  used  on  bank  angle 
with  angle-of-attack  set  to  L/D  max.  Additional  parameters  were  used  on 
angle-of-attack  when  temperature  constraints  were  considered,  as  will  be 
shown  in  Section  III. 


The  methods  considered  for  parameter  optimization  are  those  described 
in  Reference  4.  After  some  early  experimentations  with  these  parameter 
optimization  methods,  a  sequential  combination  of  two  searches  was  found 
to  produce  satisfactory  results.  This  combination  was  an  adaptive 
"creeping"  search  followed  by  a  "pattern"  search.  The  adaptive  creeping 
search  operates  by  introducing  perturbations  in  the  control  parameters 
varying  one  parameter  at  a  time  seeking  improvement  in  the  function  to  be 
optimized.  If  improvement  is  not  found  when  a  parameter  is  perturbed 
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In  one  direction  the  other  direction  Is  searched.  The  size  and  direction 
of  the  perturbations  are  determined  adaptively  by  the  search  algorithm 
as  the  search  progresses.  After  one  complete  cycle  of  the  creeping 
search,  with  perturbations  in  all  parameters  having  been  explored,  the 
pattern  search  was  used.  The  pattern  search  attempts  to  accelerate 
convergence  by  exploring  the  direction  revealed  by  the  creeping  search 
and  moves  al 1  the  parameters  simultaneously.  The  creeping  search  is  then 
repeated  and  followed  by  another  pattern  search,  etc.  Constraints  are 
handled  by  a  penalty  on  the  payoff,  forming  an  augmented  performance 
function  of  the  form: 


4>*  =  4>  +  ^  -$.|-  TT0L./2)2 


(4) 


where 

is  the  value  of  the  pay-off 
function 

Wj  is  the  weighting  matrix 

ipi  is  the  value  of  constraint  i 

iL  is  the  desired  value  of 
constraint  i 

TTOL..  is  the  tolerance 

5.  CROSS-RANGE  OPTIMIZATION 


(cross-range  used  in  this  study) 
(unit  matrix  used) 

(down-range  and/or  temperature  used) 

(down-range  2,000  -  14,000  nmi) 
(temperature  2550  -  2600°R) 

(5°  on  temperature,  50  nmi  on 
duwn-range) 


The  controls  used  to  maximize  cross-range  were  angl e-of-attack  for 
L/D  max  and  bank  angle  as  a  function  of  turn  angle.  Figure  2  shows  the 
nominal  control  used  for  the  bank  angle.  The  circles  show  the  nominal 
value  of  the  parameters  used  for  an  eight  parameter  optimization.  The 
heading  angle  points  were  selected  from  -15°  to  90°  in  15°  increments. 
The  bank  angle  at  these  points  was  then  varied  by  the  optimizer.  Linear 
interpolation  between  points  was  used.  The  first  problem  considered  was 
finding  maximum  unconstrained  cross-range.  The  down-range  constraints 
were  then  introduced  and  the  problem  re-run  seeking  maximum  cross-range 
under  a  series  of  down-range  constraints  to  develop  the  footprint.  This 
process  was  repeated  for  each  of  the  three  vehicles  studied. 
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Figure  2.  Bank  Angle  Control  Schedule 
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Typical  convergence  for  maximizing  unconstrained  cross-range  is  shown 
in  Figure  3  for  the  L/D  =  2.2  vehicle.  This  figure  shows  "steps"  in 
cross-range  at  several  places.  These  steps  occur  where  the  convergence 
process  switches  from  creep  to  pattern.  Approximately  110  iterations 
were  used  in  this  example.  Note  that  over  95%  of  the  maximum  cross-range 
was  obtained  in  30  iterations.  The  inset  to  Figure  3  shows  the  nominal 
control  schedule  used  to  start  the  optimization  and  the  final  control 
determined  by  the  optimization  process. 

This  optimized  control  and  the  resulting  controls  for  maximizing 
cross-range  under  down-range  constraints  are  shown  by  the  circles  connected 
by  the  solid  lines  in  Figure  4.  The  form  of  the  control  resulting  from 
the  optimization  suggests  that  a  linearization  might  yield  satisfactory 
results.  The  control  was  then  linearized  as  shown  by  the  dashed  lines  in 
the  figure  and  the  problem  re-run.  The  difference  in  cross-range  between 
using  the  linearized  and  the  optimal  control  was  less  than  five  percent 
in  all  cases.  A  two-parameter  linear  control  was  then  investigated.  The 
control  used  was  <J>B  =  K(a  -  af)  where  K  is  the  slope  of  bank  angle  versus 
heading,  o  is  the  instantaneous  value  of  heading,  and  is  the  final 
heading  where  the  bank  angle  reduces  to  zero.  This  control  is  illustrated 
by  Figure  5.  In  keeping  the  slope,  K,  constant  and  varying  a 
parametric  control  set  was  generated  as  shown  in  Figure  6.  As  shown  in 
the  figure,  the  maximum  bank  angle  was  arbitrarily  limited  to  75  degrees 
since  this  seemed  to  be  a  reasonable  upper  limit.  Three  different 

parametric  sets  were  run  using  K  values  of  .5,  .67,  and  1.0.  This  was 

done  for  all  three  vehicles  and  the  results  are  shown  in  Figure  7.  The 
optimum  results  are  shown  by  the  squares.  The  solid  lines  are  fairings 
through  the  outer  boundary  of  the  parametric  data.  This  figure  shows 
fairly  good  agreement  between  the  parametric  and  optimum  results  at  the 
forward  part  of  the  footprint,  but  poor  agreement  at  the  back  part.  In 
an  attempt  to  improve  the  agreement  on  the  back  part  of  the  footprint, 
the  two-parameter  control  was  changed  by:  (1)  increasing  the  bank  angle 
limit  from  75  to  85°  and  (2)  by  setting  a^.  =  0°,  which  results  in 

<pD  =  Ko.  By  varying  the  parameter  K,  a  parametric  control  set  was 

D 

generated  as  shown  in  Figure  8.  The  parameter  K  was  varied  from  0.  to  1.5. 
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Figure  5.  Two  Parameter  Control 
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Figure  6.  Two  Parameter  Bank  Angle  Control 
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The  higher  values  of  K  (K>1.0)  provide  a  bank  angle  control  which 
Initially  follows  the  85°  boundary  (Figure  8)  and  then  is  reduced  lin¬ 
early  to  zero  as  the  turning  angle  approaches  90°  (o  =  0).  This  type 
of  control  results  in  fast  initial  turning  and  short  down-range.  The 
results  are  shown  by  the  circles  In  Figure  9.  The  squares  are  again  the 
optimized  results  and  the  solid  lines  are  the  results  from  the  previous 
parametrics.  The  figure  shows  that  the  results  from  the  <{>„  =  Ka  control 

D 

compare  closely  over  the  entire  region  with  both  the  optimized  results 
at  the  short  to  mid  ranges  and  with  the  best  cross-range  achieved  by  the 
previous  parametrics  at  the  longer  ranges  for  all  three  vehicles  being 
studied.  Use  of  the  bank  control  equation  <f>B  =  K a  is  called  a  "single 
parameter"  control  as  K  is  the  only  parameter  being  varied. 

6.  FOOTPRINT  GENERATION 

The  footprint  generator  development  thus  far  has  dealt  with  maximizing 
cross-range  for  the  portion  of  the  footprint  that  has  down-ranges 
4000  nml  and  greater.  Attention  is  now  given  to  ranges  of  less  than 
4000  nmi.  Using  the  L/D  =  1.89  configuration,  the  cross-range  was 
maximized  using  the  optimization  process  where  down-range  was  constrained 
to  3000  and  then  2000  nmi.  The  resulting  optimal  bank  angle  controls 
are  shown  in  Figure  10  for  these  ranges  plus  the  4000  nmi  range  developed 
previously.  It  was  found  that  these  results  could  be  closely  approximated 
using  the  linear  control  4>g  =  Kc  if  the  maximum  bank  angle  limit  were 
removed,  i.e.,  allowing  the  vehicle  to  bank  greater  than  85  degrees. 

The  results  are  shown  in  Figure  11.  Note  that  all  trajectory  data  points 
were  run  at  (L/D)  max  except  for  the  zero  lift  data  point.  Complete 
footprints  for  the  three  configurations  are  shown  in  Figure  12.  This 
figure  shows  good  agreement  between  the  optimal  results  shown  by  the 
hexagonal  points  and  the  single  parameter  control  solutions. 

From  Figure  12  each  footprint  was  generated  with  13  points,  each 
point  requiring  a  trajectory.  In  contrast,  each  of  the  optimization 
points  shown  in  the  figure  represents  about  110  trajectory  iterations 
before  converging  to  a  solution.  Forming  part  of  a  footprint  using  three 
optimized  points  would  require  about  110  X  3  =  330  trajectories. 
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Figure  9.  Optimal  and  Parametric  Comparisons 
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Figure  10.  Optimal  Bank  Angle  Schedules 


Figure  11.  Comparison  of  Optimal  and  Parametric  Data 
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Use  of  the  linear  control  to  produce  a  complete  footprint  hence  requires 
less  computer  resources  by  a  factor  of  330/13  =  25.  Converting  this  to 
computer  solution  time  using  an  average  of  30  seconds  computer  time  per 
trajectory  on  a  CDC  CYBER-175  yields  a  decrease  in  time  from  9900  seconds 
for  the  three  optimized  points  to  only  390  seconds  for  the  complete 
footprint  using  the  linear  control. 

Summarizing  then,  developing  a  footprint  by  the  single  parameter 
control  procedure  was  found  to  closely  approximate  optimization  results 
for  the  vehicles  studied.  The  procedure  outlined  below  is  for  an  initial 
equatorial  orbit  heading  east.  For  other  inclinations,  z  (remaining 
heading  change)  must  be  substituted  for  a  in  the  bank  angle  control,  where 
Z  is  the  angle  between  the  component  of  the  velocity  vector  parallel  to 
the  earth  and  the  cross-range  axis.  Single  parameter  control  procedure: 

a.  The  Initial  heading  and  position; are  set.  For  the  equatorial  orbit 
heading  east,  o  s  90  degrees  and  the  position  is  fixed  on  the  equator. 

V. 

b.  The  independent  variable  for  the  bank  angle  control  is  heading, 
o  (or  In  general  z)  ■  Two  points  are  required,  one  at  the  beginning  (aQ) 
and  one  at  the  end  (a^)  of  the  turn  (for  initial  equatlonal  orbits  0^=0), 
Linear  interpolation  between  these  points  forms  the  bank  control. 

c.  The  initial  bank  angle  (at  oQ)  is  varied  parametrically  in 
15-degree  increments  from  0  to  150  degrees  to  form  the  footprint.  The 
final  heading  angle  a ^  is  0;  hence  the  term  "single  parameter  control"  as 
only  the  bank  angle  at  aQ  is  being  varied  to  form  the  linear  control. 

d.  One  more  point  (minimum  down-range)  Is  needed  to  complete  the 
footprint.  In  this  study  the  minimum  down-range  point  was  run  at  zero 
lift  (pull  downs  were  not  considered). 

7.  SOME  COMPARISONS  AND  PERTURBATIONS 

The  remaining  work  to  be  discussed  Involves  the  comparison  of  various 
entry  conditions  and  control  schedules.  Figure  13  shows  how  the  footprint 
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Figure  13.  Effects  of  Bank  Initiation  on  Footprint 
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changes  due  to  initiating  bank  at  the  beginning  of  the  sensible  atmosphere, 
295K  ft,  and  at  pull-out,  y  =  0.  Both  footprints  were  generated  using 
the  single  parameter  control  schedule  and  the  L/D  =  2.2  vehicle. 

A  bank  angle  control  schedule  much  used  in  the  past  for  footprint 
generation  was  to  hold  a  constant  bank  angle  until  the  heading  changed 
by  90  degrees  and  then  reduce  the  roll  to  zero.  By  varying  the  amount 
of  bank,  a  footprint  is  generated.  Figure  14  compares  a  footprint 
using  this  type  of  control  with  the  single  parameter  control.  For  this 
case  the  L/D  =  1.89  vehicle  was  used. 

The  effect  of  initial  path  angle  on  the  footprint  area  available  is 
shown  in  Figure  15  for  the  L/D  =  1.89  vehicle.  The  single  parameter 
control  was  used  to  generate  this  plot.  The  corresponding  pull -cut 
altitudes,  y  =  0,  are  shown  in  Figure  16.  It  should  be  noted  that  for 
a  given  set  of  entry  conditions  the  initial  pull-out  altitude  is  a 
function  of  the  parameter  W/(CLS)  rather  than  L/D.  The  Figure  16  data 
are  for  entry  at  L/D  max,  and  the  altitudes  shown  are  approximately  the 
same  for  the  L/D  =  2.2  and  L/D  =  2.6  vehicles.  The  pull-out  altitude  can 
vary  appreciably  depending  upon  the  angle-of-attack  schedule  used  during 

'V 

entry.  The  lower  the  value  of  W/(C^S)  the  higher  the  initial  pull-out 
altitude  will  be.  This  illustrated  by  Figure  17  which  shows  entry 
trajectories  flown  both  at  max  and  at  (L/D)  max  to  the  pull-out  point, 
y  =  0.  This  figure  shows  that  the  max  entry,  having  a  lower  value  of 
W/(CLS),  results  in  pull-out  at  an  appreciably  higher  altitude  than  the 
(L/D)  max  trajectory.  Range  is  also  increased  for  the  max  entry  for 
the  case  where  altitude  is  held  constant  after  pull-out. 

Trajectories  that  are  flown  at  L/D  max  or  at  a  constant  angle-of-attack 
typically  oscillate  in  altitude.  These  oscillations  can  be  damped  using 
bank  angle  and/or  angle-of-attack  control.  To  maintain  leve"  flight  after 
pull-out,  the  lift  is  reduced  such  that  the  lift  plus  the  centrifugal 
force  is  equal  to  the  weight.  For  straight  down-range  (no  bank)  the 
angle-of-attack  is  reduced  to  the  proper  value  at  pull-out  and  is  increased 
as  the  speed  decreases  due  to  aerodynamic  drag.  The  angle-o^-attack  for 
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Figure  14.  Comparisons  of  Single  Parameter  and 
Constant  Bank  Angle  Controls 
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L/D  max  is  flown  once  the  L/D  max  equilibrium  glide  trajectory  is  reached. 
In  the  equilibrium  glide  trajectory  the  vertical  forces  (lift  plus 
centrifugal  force)  nearly  equal  the  gravitational  forces  and  hence  the 
vehicle  descends  without  altitude  oscillations.  One  disadvantage  of  this 
type  of  trajectory  is  a  great  reduction  in  range  since  only  a  small  part 
of  the  profile  is  flown  at  L/D  max.  Other  variations  for  smoothed 
trajectories  would  be  to  raise  the  pull-out  altitude  by  flying  at  CL  max 
until  pull-out  and/or  flying  a  very  small  positive  or  negative  path  angle 
after  pull-out.  The  effects  of  flying  small  controlled  flight  path  angles 
after  pull-out  are  shown  in  Figure  18.  This  figure  also  shows  the  (L/D) 
max  trajectory  for  comparison. 

Cross-range  can  be  obtained  either  at  constant  altitude  or  at  small 
flight  path  angles  afer  pull-out  by  using  bank  angle  control  and  (L/D) 
max  for  the  turn.  By  rolling  out  to  wings  level  and  switching  to 
angle-of-attack  control  once  a  given  heading  change  is  reached,  a 
footprint  can  be  generated,  as  shown  by  Figure  19.  The  maximum  amount  of 
turn  available  is  to  maintain  the  turn  using  bank  control  until  the  wings 
are  level.  This  is  shown  by  the  dashed  lines  in  Figure  19.  The  altitude 
profiles  for  the  down-range  points  with  no  turn  in  Figure  19  are  the 
same  as  those  shown  in  Figure  18  for  corresponding  values  of  the  flight 
path  angle  control,  y  . 

The  single  parameter  control  and  the  constant  altitude  control  have 
been  used  to  generate  a  footprint  data  base.  These  data  are  shown  in 
Appendix  C. 
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Figure  18.  Comparison  of  Smoothed  and  (L/D)  Max  Trajectories 


Figure  19.  Lateral  Range  Using  Bank  Angle  Control 
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SECTION  III 

IMPACT  OF  HEATING  CONSIDERATIONS 

In  an  effort  to  assess  the  Impact  of  aerodynamic  heating  constraints 
on  entry  vehicle  maneuvering  capability,  the  low  (L/D)  configuration  was 
selected  for  study.  The  temperatures  discussed  in  this  section  were 
computed  at  a  point  approximately  ten  feet  aft  of  the  nose  on  the  lower 
surface  centerline  of  the  vehicle  using  a  one-dimensional  conduction 
analysis  on  a  simplified  structure. 

1 .  APPROACH  AND  PROCEDURES 

The  approach  used  was  to  calculate  the  heat  input  to  the  surface 
through  the  use  of  a  "reference”  heat  transfer  coefficient  which  accounted 
for  the  location  of  the  vehicle  on  the  entry  trajectory  (altitude  and 
velocity)  and  a  ratio  of  local -to-reference  heating  coefficient  to  account 
for  vehicle  attitude  and  point  of  interest  on  the  vehicle.  This  approach 
is  similar  to  that  used  in  Reference  7.  The  following  heating  relationships 
were  used  to  compute  the  radiation  equilibrium  wall  temperature. 


(h/hF&R)  =  f(a) 

(5) 

q  =  hp&R  (h/hp&R)  (Haw  -  Hw) 

(6) 

T  =  Cq/ (4.75  X  10‘13  x  e)]  ,25 

(7) 

where  the  ratio  (h/hpjR)  is  obtained  from  theory  or  wind  tunnel  data 
correlations  and  ( hp&R)  is  calculated  using  the  Fay  and  Riddell  theory. 
Values  of  (h/hp&R)  with  vehicle  attitude  used  in  this  study  are  shown  in 
Table  1. 

To  maintain  continuity,  the  same  parameter  optimizat'on  techniques 
used  previously,  namely  "creep"  and  "pattern",  were  used  for  the  heating 
analysis.  To  provlce  the  capability  for  constraining  the  equilibrium 
wall  temperature  to  specified  values,  the  same  penalty  function  approach 
used  previously  for  the  down-range  constraint  was  employed. 
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TABLE  1 

HEATING  VARIATION  WITH  VEHICLE  ATTITUDE 


(h/hF&R) 

a 

.17 

0 

.175 

20 

.20 

25 

.23 

30 

.23 

40 

.26 

45 

.30 

50 

The  control  philosophies  employed  to  generate  the  unconstrained  and 
temperature  constrained  reentry  footprints  are  outlined  in  Figure  20. 

For  the  unconstrained  footprint  (no  heating  constraints),  problems  were 
formulated  whereby  the  vehicle  flew  an  angle-of-attack  schedule  for 
maximum  lift-to-drag  ratio  over  the  entire  reentry  trajectory.  A  oank 
angle  of  zero  degrees  was  maintained  from  entry  to  the  initial  flare 
point  (y-0  degrees).  From  the  initial  flare  point  to  trajectory 
termination  (H=50,000  feet),  the  program  modulated  the  bank  angle  to 
maximize  the  cross-range.  For  the  temperature  constrained  problems,  both 
of  the  control  parameters  were  modulated  over  the  entire  trajectory  with 
the  exception  that  during  the  initial  plunge  into  the  atmosphere  the  bank 
angle  was  set  to  zero  degrees.  Using  this  steering  philosophy  ensured 
that  the  heating  criteria  would  not  be  violated  at  any  point  along  the 
entry  flight  profile.  Listed  in  Table  2  are  the  initial  entry  conditions 
used  for  the  aerodynamic  heating  analysis. 
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Figure  20.  Control  Philosophies  Used  for  Reentry  Trajectories 
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TABLE  2 

INITIAL  ENTRY  CONDITIONS 
hQ  =  364,800  feet 

V.  =  26,031  feet  per  second  (Inertial) 

o 

VQ  =  24,480  feet  per  second  (Relative) 

Y0  *  -2.0  degrees 

aQ  =  90  degrees 

*0=00=0  degrees 
(W/S)  =  70  #/ft2 

All  of  the  trajectories  were  terminated  at  50,000  feet  altitude. 

To  investigate  the  Impact  of  aerodynamic  heating  constraints  on 
reentry  vehicle  maneuverability.  It  was  desired  to  generate  a  footprint 
with  temperature  constraints  to  compare  to  the  unconstrained  footprint. 
Shown  in  Figure  21  is  the  unconstrained  reentry  footprint  for  the  low 
llft-to-drag  ratio  configuration.  Except  for  the  maximum  down-range 
point,  all  of  the  points  on  the  footprint  were  obtained  by  maximizing 
the  lateral  range  at  various  down-range  distances.  Although  the 
radiation  equilibrium  temperature  (RET)  was  not  constrained  It  was 
computed  on  all  of  the  reentry  trajectories.  The  maximum  value  of  cross¬ 
range  (1964  nmi)  was  attained  at  approximately  6000  nmi  down-range.  For 
this  case  the  maximum  radiation  equilibrium  temperature  was  2670°R. 

2.  RESULTS 

The  Initial  problem  Investigated,  with  temperature  constraints,  was 
the  maximum  cross-range  case.  2670°R  was  considered  to  be  the  upper* 
limit  since  that  value  was  obtained  for  the  unconstrained  case.  Figure  22 
Illustrates  the  maximum  cross-range  sensitivity  to  changes  In  the  maximum 
allowable  surface  temperature.  As  noted  In  Figure  22  two  cases  were  run 
with  temperature  constraints  of  2550  and  2600°R.  The  curve  shows  that 
decreasing  the  temperature  from  2670°R  to  2550°R  results  In  a  9.4  percent 
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Figure  21.  Reentry  Footprint  (No  Temperature  Constraints) 
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decrease  In  cross-range.  The  time  histories  of  the  surface  temperature, 
angle-of-attack,  and  bank  angle  are  shown  in  Figures  23  through  28  for 
the  three  cases.  Analysis  of  the  results  indicated  that,  as  the 
temperature  constraint  used  In  this  analysis  becomes  more  stringent,  the 
vehicle  is  forced  to  enter  at  a  lift-to-drag  ratio  considerably  less  than 
maximum,  approaching  (C^)  max  for  the  2550°R  constraint  case.  Also,  the 
shape  of  the  entry  flight  profile  is  highly  Influenced  by  the  temperature 
constraint  in  that  the  initial  pull-out  altitude  occurred  approximately 
on  the  constraint  boundary.  To  show  the  impact  of  temperature  constraints 
on  the  complete  footprint,  the  most  stringent  temperature  limit,  namely 
2550°R  was  used  in  conjunction  with  down-range  constraints  to  generate 
the  reduced  entry  footprint  shown  In  Figure  29.  Inspection  of  these 
entry  footprints  shows  the  large  amount  of  geographic  coverage 
that  is  lost  as  the  result  of  the  temperature  constraints.  The  reason 
for  this  reduction  in  maneuverability  is  that  the  vehicle  reentering 
the  earth's  atmosphere,  flying  a  temperature  constrained  flight  profile, 
is  forced  to  reenter  at  a  higher  altitude  than  the  unconstrained  case, 
therefore  a  lower  dynamic  pressure  and  lift/drag  ratio. 

It  is  interesting  to  note  that  although  the  maximum  cross-range  of 
the  footprint  was  reduced  only  184  nml  out  of  1964  when  the  maximum 
temperature  was  reduced  from  the  unconstrained  condition  to  2550°R, 
significant  amounts  of  down-range  were  lost.  On  the  forward  part  of  the 
footprint  the  maximum  down-range  was  reduced  from  14500  to  11700  nml 
and  on  the  rear  portion  of  the  footprint  the  minimum  range  Increased 
from  around  1000  to  2000  nml.  The  footprint  collapses  towards  the 
center  but  not  as  much  as  one  might  expect  considering  that  the  limit  of 
2550°R  is  about  the  lowest  that  can  be  achieved  on  the  initial  plunge 
from  the  initial  entry  conditions  outlined  in  Table  2. 

3.  SUMMARY  AND  COMMENTS 

In  conclusion,  the  footprint  generated  by  Imposing  a  temperature 
constraint  of  2550°R  compared  to  the  unconstrained  footprint  of  Section  II 
loses  9.4%  of  the  maximum  cross-range  capability,  19%  of  the  maximum 
down-range,  and  the  minimum  down-range  is  doubled  from  1000  to  2000  nmi . 
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Figure  24.  Unconstrained  Case,  Angle-of-Attack  and  Bank  History 


AFWAL-TR-82-3036 


l  u  -  1.89  y  =  -2 

.  e 


3000  r 


2500 


2000 


1500 


1000 


500 


l  £°co^»c^ooc 


—/?rr-  ttsco - TEMPERATURE  BOUNDARY 


o 

o 

o 

o 

o 

o 

o 


Oo 


(2600°R) 


o 

c. 


o 

— o - *- 

0  400 


— i _ i _ i _ i _ i 

800  1200  1600  2000  2400 

TIME  -  SECONDS 


2800 


Figure  25.  Allowable  Temperature  2600°R,  Surface  Temperature  History 
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Figure  26.  Allowable  Temperature  2600® R,  Angle-of-Attack  and  Bank 
story 
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Figure  27.  Allowable  Temperature  2550°R,  Surface  Temperature  History 
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Figure  28.  Allowable  Temperature  2550°R,  Angle-of-Attack  and  Bank 
Histories 


Figure  29.  Sensitivity  of  Temperature  Constraint  on  Reentry  Footprint 
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Due  to  the  substantial  amount  of  computer  time  Involved  to  run  optimal 
solutions,  constrained  footprints  were  not  generated  for  other  initial 
conditions.  The  heating  method  used  In  the  trajectory  model  has  provision 
for  simultaneously  constraining  temperature  on  more  than  one  point  on 
the  vehicle.  However,  the  lack  of  data  for  the  heating  model  on  the  other 
points  of  the  low  L/D  vehicle  or  the  other  vehicles  prevented  running 
them  with  temperature  constraints. 


47 


AFWAL-TR-82-3036 


SECTION  IV 

APPLICATION  AND  DISPLAY 

The  work  shown  to  this  point  has  been  to  develop  a  procedure  for 
rapid,  accurate,  footprint  calculations  under  the  various  flight 
constraints.  However,  to  determine  the  utility  of  a  vehicle  concept  in 
performing  a  given  function,  the  footprint  data  must  be  related  and 
applied  to  that  function.  For  such  concepts  as  an  orbital  based 
maneuvering  entry  vehicle,  questions  of  Interest  might  include  the 
following:  (1)  What  Is  the  landing  coverage  from  a  specified  deorbit 
point?  (2)  What  are  the  effects  of  node  rotations  and  inclination  change? 
(3)  What  portions  of  which  orbits  allow  landing  in  a  set  of  designated 
landing  sites?  A  similar  set  of  questions  could  be  phrased  for  a  boost 
glide  missile  system.  The  procedure  developed  for  the  display  and  analysis 
portion  of  the  footprint  generator  was  designed  as  a  tool  to  find 
answers  to  those  questions. 

Options  are  provided  to  display  the  footprint  on  two  types  of  world 
maps  and  calculate  a  modified  footprint  to  allow  the  user  to  change  the 
node  location,  orbital  Inclination,  or  increase  down-range  for  stepping 
the  deorbit  point  around  an  orbit.  An  option  Is  also  provided  to  allow 
input  of  only  half  of  a  symmetrical  footprint  and  computer  generate  the 
other  half.  Required  Inputs  to  the  procedure  are  longitude,  latitude, 
and  heading  angle  of  a  reference  point  on  the  Initial  orbit  and  the 
longitude,  latitude  points  of  the  previously  calculated  footprint.  If  a 
modified  footprint  Is  to  be  calculated,  the  procedures  shown  below  on 
"calculating  a  new  footprint"  are  used  first.  If  display  only  of  the 
generated  footprint  Is  desired  this  section  Is  skipped. 

1.  CALCULATING  A  NEW  FOOTPRINT 

The  first  step  Is  to  take  the  longitude-latitude  points  from  the 
footprint  boundary  calculated  by  the  footprint  generator  and  convert  this 
data  into  cross-range  and  down-range  from  a  reference  such  as  the  deorbit 
condition  (longitude,  latitude  of  the  deorbit  point  plus  the  heading 
angle  of  the  orbit  at  that  point  to  define  the  initial  great  circle  for 
down-range  calculation). 
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The  geometry  of  the  problem  is  shown  In  Figure  30  where  point  L  is 
the  initial  point,  arc  LP  is  the  initial  great  circle,  and  point  F  is  a 
point  on  the  footprint  boundary.  The  cross-range  of  point  F  is  defined 
as  the  perpendicular  distance  from  the  point  to  the  initial  great  circle, 
shown  in  the  figure  as  the  arc  PF.  The  down-range  is  then  the  distance 
along  the  initial  great  circle  from  the  initial  point  to  the  point  at 
which  the  cross-range  Is  measured  (arc  LP). 

From  spherical  trigonometry,  the  great  circle  range  LP  to  the  point  F 
Is  computed 

LF  =  cos"^[sin  4>l  sin  <|>Lq  +  cos  ^  cos  <^q  cos  (0^  -  0L^ ) ]  (8 
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The  heading  c  of  the  great  circle  from  the  Initial  point  to  the 
point  F  Is  computed  from  the  spherical  triangle  LNF: 

tan  c  =  sin  (0L  -  0Lg)cos  <J>l  cos  <j>LQ 
sin  <{>L  -  sin  <j>LQ  cos  (LF) 


and  the  triangle  LPF  Is  solved  for  the  down-range  Xp,  and  cross-range  Yj. 

Xp  =  R'  cos  _1  [cos  (LF)/cos(s1n”^  (sin(LF)  sin  £))]  (1C 

Yp  =  R1  sin  [sin  (LF)  sin  5]  (11 


where:  £  =  £  -  aQ 

and  R*  Is  the  average  radius  of  the  earth 

The  next  step  is  to  modify,  through  Input,  the  node  location  (by 
changing  G^p  of  the  Initial  point),  the  Inclination  angle  of  the  Initial 
orbit  (by  changing  Op),  or  the  down-range  (Xp  which  was  just  calculated 
for  stepping  around  the  orbit).  These  parameters  may  be  changed  singly 
or  In  combination. 

Next,  this  revised  condition  Is  then  converted  back  to  a  new 
longitude,  latitude  point  for  plotting  and  analysis.  However,  what  was 
a  straightforward  calculation  going  to  cross-range  and  down-range  now 
requires  an  interatlon  converting  back  to  longitude  and  latitude. 
Starting  from  the  new  down-range  (XpM)  and  heading  at  the  initial  point. 
(opM)  holding  the  cross-range  (arc  PF)  constant  and  rearranging 
Equations  3  and  4: 

LFm  =  cos-1  [cos  (PF)  cos  (XpM/R*)]  (12 

where  LFM  is  the  new  value  of  arc  LF 

£  M  *  sin-1  [sin  (PF)/s1n  (LF^)]  (13 

ind  CM  ’  *  °0M  <14 


where  LFM,  £M,  and  cM  are  the  new  values  of  LF,  £  and  c. 
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Now  a  guess  at  a  new  latitude  of  the  footprint  point  Is  made  (4^), 
and  (6^  -  9lq)m  Is  computed  from  Equation  8.  Putting  this  In  Equation  9 
yields  a  5  which  Is  compared  with  the  result  obtained  from  Equation  14. 
If  the  values  of  c  and  are  not  equal,  a  new  4^  Is  chosen  and  the 
process  repeated  until  It  converges. 


The  new  longitude  Is  then  computed  by: 

0LM  =  (0L  "  W  M  +  0  LOM 


(15) 


where  the  subscript  M  denotes  tho  new  value  of  the  variables.  This 
procedure  Is  repeated  until  all  points  of  the  footprint  have  been 
processed.  At  this  point  the  new  footprint  Is  available  for  plotting. 


2.  PLOTTING  A  FOOTPRINT 

Display  of  the  footprint  Is  accomplished  by  plotting  the  longitude 
and  latitude  points  on  a  map  from  the  world  mapping  capability  available 
through  the  DISSPLA  Integrated  software  system  and  plotting  language 
(Reference  12).  From  the  15  map  projections  available  In  this  system 
two  were  chosen  to  display  the  results,  the  cylindrical  equidistant  and 
the  orthographic  projection.  All  or  any  portion  of  the  world  can  be 
projected.  The  cylindrical  equidistant  shows  the  coordinates  as  they 
are  on  a  flat  map  and  In  the  orthographic  projection  one  views  the  globe 
from  far  away-projectlon  pole  at  Infinity. 


The  plotting  program  is  structured  to  make  both  these  maps  for  each 
condition  desired  to  be  plotted.  As  only  half  the  earth  Is  visible  on 
the  orthographic  projection,  the  viewing  angle  Is  varied  so  that  all  the 
footprint  (or  most  for  large  footprints)  can  be  seen  unless  the  user 
Instructs  otherwise.  Some  sample  results  as  presented  In  the  fo' lowing 
paragraphs  were  copied  from  a  Tektronix  terminal. 

Figure  31  shows  the  footprint  for  the  L/D  =  1.89  vehicle  when  operated 
at  a  high  angle-of-attack  to  reduce  heating  loads.  The  footprint  was 
run  with  the  starting  point  taken  as  the  entry  conditions  of  Table  2  and 
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Figure  31.  Sample  Plots  L/D  =  1.39  Reduced  Footprint 
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heading  due  east  at  0°  latitude  and  longitude.  Note  the  starting  point 
did  not  lie  on  the  footprint.  The  positive  latitude  half  of  the  footprint 
was  Input  and  the  other  half  was  calculated  by  symmetry  in  the  display 
program.  The  left  side  of  the  figure  shows  plotting  on  the  cylindrical 
equidistant  projection.  At  the  top  in  Figure  31  (a)  is  the  footprint 
"as  run"  at  a  heading  angle  of  90°  (inclination  0°),  the  starting  point  at 
0°  latitude  and  longitude  is  shown  by  the  filled  circle.  In  Figure  31  (b) 
the  initial  heading  has  been  changed  to  30®  (inclination  60°),  and  at 
the  bottom  in  Figure  31(c)  a  heading  of  -90®  is  shown.  The  right  side 
of  the  figure  shows  inclination  0°  in  Figure  31  (d)  and  inclination  60® 
in  Figure  31(e)  on  the  orthographic  projection. 

Figure  32  shows  similar  plots  for  the  larger  footprint,  this  for  the 
L/D  =  2.6  vehicle.  The  starting  point  was  again  the  conditions  from 
Table  2  and  heading  east  at  0®  longitude  and  latitude.  The  left  side  of 
the  figure  on  the  cylindrical  equidistant  projection  shows  inclination 
of  0°  in  Figure  32  (a)  and  30°  in  32  (b).  Note,  that  the  points  are  not 
connected  when  the  lines  would  cross  longitude  180°,  The  right  side 
shows  the  same  plots  on  the  orthographic  projection.  The  area  of  the 
footprint  is  so  large  it  won't  fit  on  half  the  globe  so  only  a  portion 
can  be  shown  on  the  orthographic  projection.  The  area  of  the  map  not 
covered  by  the  footprint  is  shaded  for  ease  of  viewing.  Figure  32(c) 
shows  the  inclination  0®  and  32(d)  inclination  30°. 

3.  EARTH  ROTATIONAL  EFFECTS 

An  approximation  to  earth  rotational  effects  has  been  included  as  an 
option  in  the  computer  model.  To  correctly  include  rotation  it  would  be 
necessary  to  re-run  trajectories  for  each  case  being  considered.  However, 
if  it  is  assumed  that  the  basic  trajectory  is  unchanged  by  a  rotating 
earth  and  the  predominate  effect  is  that  of  the  earth  rotating  below  the 
trajectory  a  first  order  approximation  of  the  shift  due  to  the  rotation 
can  be  obtained.  This  is  a  shift  in  the  longitude  of  the  end  point,  with 
the  latitude  unchanged,  and  is  computed  by  6LRqT  =  6LNR0T~a)pt’’  where 
0LROT  *s  t*ie  new  ^on9ltuc*e  the  encl  P°*nt  1n  degrees,  9|_NRqt  is  the 
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longitude  of  the  nonrotating  and  point  Is  earth  angular  rotation 
rate  =  .004178  deg/sec,  and  t  Is  the  flight  time  of  the  trajectory.  The 
flight  time  (t)  Is  obtained  from  the  nonrotating  trajectory  run';  which 
developed  the  footprint. 

Typical  flight  times  range  from  2000  to  4000  seconds  and  are  generally 
shorter  on  the  back  portion  of  the  footprint  than  they  are  on  the 
forward  portion.  This  then  produces  a  longitude  shift  of  8.4  to  16.7 
degrees.  Note  that  the  longitude  shift  is  toward  the  west,  that  is  toward 
more  negative  values.  This  has  the  effect  of  stretching  the  footprint 
when  the  initial  conditions  are  heading  west  on  the  equator,  compacting 
it  when  heading  east,  and  distorting  it  at  other  inclinations  as  shown 
in  Figure  33.  The  starting  point  in  Figure  33  taken  as  the  ent*y 
conditions  of  Table  2,  was  at  0°  longitude,  0°  latitude  and  is  ^ut  of 
view  in  the  diagram.  This  point  corresponds  to  the  time  =  0  condition 
and  would  not  be  moved  by  earth  rotation. 


Figure  33.  Shift  Due  to  Earth  Rotation 
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4.  OTHER  PLOTTING  CONSIDERATIONS 

This  system  as  described  represents  an  attempt  to  automate  a  procedure 
for  display  and  manipulation  of  footprints  requiring  little  user  knowledge 
of  the  DISSPLA  system  and  a  minimum  of  Input  data.  Many  other  different 
things  may  be  done  to  tailor  a  program  to  produce  the  most  useful  plots 
for  a  particular  application.  Among  these  would  be:  other  types  of 
projections,  mapping  of  only  a  portion  of  the  earth,  labeling  and  titling 
plots,  shading  of  footprints.  Insertion  of  vehicle  pictures,  display  of 
more  than  one  footprint  on  the  same  map,  or  structuring  of  the  analysis 
program  to  generate  and  plot  mission  oriented  features.  Many  of  these 
can  be  done  with  very  little  effort. 

The  footprint  generator  developed  in  this  report  shows  potential  as 
a  powerful  tool  for  rapid  calculation  of  entry  vehicle  maneuvering 
capability.  The  application  and  display  work  developed  here  provides  a 
first  step  towards  rapidly  analyzing  vehicle  application  to  a  required 
mission  and  producing  high  quality  plots. 


56 


AFVIAL-TR-82-3036 


SECTION  V 
CONCLUSIONS 

1.  Parameter  optimization  methods  were  applied  to  maximize  lateral  range 
with  and  without  downrange  constraints  to  produce  performance  footprints. 
Good  convergence  was  obtained  from  using  the  methods  "creeper"  and 
"pattern"  on  the  bank  angler  control  schedule,  angle-of-attack  was  set  to 
L/D  max. 

2.  A  simplified  bank  control  procedure  was  developed  which  closely 
approximated  the  optimal  footprint  results  where  no  heating  constraints 
were  Imposed.  The  procedure  Is  outlined  in  Section  II.  The  method  Is 
easy  to  use  and  Is  efficient  requiring  less  computer  resources  than  the 
optimization  procedure  by  a  factor  of  about  25  to  1 . 

3.  No  simple  method  was  found  to  handle  heating  constraints!.  However, the 
parameter  optimization  procedure  was  able  to  produce  solutions  for  the 
single  heating  point  considered.  Footprint  shrinkage  when  applying  the 
temperature  constraint  was  most  noticeable  at  the  front  and  back  portions 
with  less  shrinkage  In  maximum  cross-range.  For  the  most  stringent 
temperature  constraint  considered  the  results  were  a  19%  loss  in  maximum 
down-range,  a  doubling  of  the  minimum  down-range  and  a  9.4%  loss  In 
maximum  cross-range. 

4.  An  automated  procedure  was  developed  to  display  and  move  the  footprint 
on  world  maps.  This  technique  shows  promise  in  relating  results  to 
functional  performance  and  can  be  used  separately  or  in  conjunction  with 
the  footprint  data. 
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APPENDIX  A 

COORDINATE  SYSTEM  TRANSFORMATIONS 

This  appendix  develops  the  transformations  necessary  to  relate  the 
coordinate  systems  employed. 


Local -Geocentric 
Horizon  Coordinates 


X,  Y,  Z 


Earth-referenced  coordinates 
Local -geocentric -horizon  coordinates 


To  describe  body  motion  relative  to  the  planet,  a  local-geocentric- 

To  locate  the  X  ,  Y 


horizon  coordinate  system  was  employed. 

'e- 


axes 

g  g  g 

with  respect  to  the  Xg,  Yg,  Zg  axes,  rotation  was  first  about  Zp  by  an 
angle  (180°  +  C)  and  then  about  Y^  through  the  angle  (90°  -  <f>^)  where 
is  vehicle  latitude  and  the  angle  C  is  related  to  the  vehicle  position 
by  the  expression: 


C  =  tan 


-1 


(Ye/Xe) 


(Al) 
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The  complete  transformation  reduced  to  a  single  matrix  is  given  by: 


-s1n(<}>L)cos{c) 

sln(c) 


-COS  (  <t>L )  cos  ( c ) 


-s1n(<j>L)s1n(c) 

-cos(c) 


-cos  (<J>L)sin{c) 


-cos(d>L) 


sin(<pL) 


Aerodynamic  forces  are  conveniently  summed  in  a  wind  axis  coordinate  system 
(XA,  Y^,  ZA).  Since  the  equations  of  motion  are  solved  In  (Xe>  Yfi,  Zg)  coordinates 
the  wind  axis  components  of  force  must  then  be  resolved  Into  this  basic  system. 

The  transformation  from  local -geocentric-horizon  coordinates  to  wind  axes  is  given  by: 


cos(ya)cos (aA) 
-sinfa^cosf^) 


cos(yA)sin(aA) 
cos(o.)cos(4>  ) 


-sin(yA) 


+s1n(yA)cos(aA)s1n(<t>B)  +sin(aA)s1n(aA)sin($  ) 


sin(oA)sin(<t>B) 

+sin(yA)cos(oA)cos(<f>B) 


-cos(aA)sin{<|>B) 

+sin(yA)sin(aA)cos(<t>B) 


cos(yJsin(<t>  )  |  Y 

M  B  i  y 


cos(aA)cos(^B)  Zg 


which  defines  a  direction  cosine  set 


the  resolution  from  wind  axes  to  local-geocentric  then  become 
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The  local -geocentric  components  are  then  resolved  Into  the  Xg,  Y  ,  Zg 
system  through  a  transposition  of  the  1,j,k  matrix,  as  follows: 


The  forces  acting  on  the  vehicle  are  the  sum  of  the  aerodynamic  and 
gravity  forces. 

EFX  ■  D  ♦  m9x 

EFy  •  Y  -  0  (AS) 

IFZ  =  L  +  mgz 

The  side  force,  Y,  is  zero  since  only  coordinated  turns  with  zero 
side  slip  were  considered.  The  drag  and  lift  terms  are: 


D  =  Cp  g  S 
L  =  CL  q  S 


(A9) 


Where  the  dynamic  pressure,  q.  Is  a  function  of  velocity  and  altitude 
and  the  lift  and  drag  coefficients  are  functions  of  Mach,  angle-of-attack, 
and  altitude.  Thrust  was  not  used  In  this  analysis. 
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where 

P2  ■  1  -  3  sin2  <j>L 

3 

Pj  a  3  sin  <f>^  -  5  sin 

P4  =  3  -  30  sin2  <f>^  +  35  sin*  4^ 

P&  *  sin  <J>L  cos  <j»L 

2 

Pg  *  cos  <j>L  (1  -  5  sin  4>l) 

2 

?7  =  sin  4>l  cos  4>l  (-3  +  7  sin  <|>L) 

The  constants  used  for  the  planet  earth  are  as  follows: 

y  =  1.407698  x  101 6  ft3/sec2 
R  =  20,925,631  ft 
J  =  1623.41  x  10'6 
H  *  6.04  x  10"6 
K  ■  6.37  x  10“6 

The  polar  radius  used  In  the  oblate  earth  calculations  for  latitude 
and  altitude  Is  Rp  =  20,885,965  ft. 
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APPENDIX  B 

AERODYNAMIC  CHARACTERISTICS 

This  appendix  presents  the  aerodynamic  coefficient  data  used  for 
trajectory  calculations.  Sketches  of  the  configurations  are  shown  In 
Figure  la  through  lc  In  the  main  part  of  this  report  and  In  Figures  C-l 
through  C-3  In  Appendix  C.  'These  data  were  obtained  from  References  8,  9, 
10,  13,  and  14. 

Tables  B-l  and  B-2  present  the  aerodynamics  for  the  vehicle  shown  in 
Figure  la.  Table  8-3  for  the  vehicle  of  lb,  and  Tables  B-4  and  R-5  for 
the  vehicle  of  lc.  Tables  B-6  and  B-7  show  the  aerodynamics  fo~  the 
vehicle  shown  in  Figure  C-l  and  Table  B-8  for  the  vehicle  of  C-3. 

The  configuration  shown  In  Figure  C-2  is  the  same  as  the  configuration 
shown  In  Figure  lc  in  the  main  part  of  this  report  and  its  aerodynamics 
were  given  In  Tables  B-4  and  B-5. 


(L/D)  =  1.89  CONFIGURATION  DRAG  COEFFICIENT  DATA 
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APPENDIX  C 
REENTRY  APPLICATIONS 


C-l .  INTRODUCTION 

Two  different  sets  of  data  are  presented  to  show  the  maneuvering 
capability  available  for  lifting  entry  when  two  different  types  of 
controls  are  used.  The  first  shown  as  data  set  1  below,  uses  the  single 
parameter  control  <|>B  =  Ka  discussed  in  the  main  part  of  this  report.  The 
second  shown  as  data  set  2  uses  a  combination  of  bank  and  angle-of-attack 
control  to  hold  altitude  constant  after  pull-out. 

Three  vehicles  were  used  for  this  analysis:  (1)  a  conventional  circular 
missile  having  "low"  (L/D)  characteristics,  (2)  a  reentry  configuration 
having  "medium"  (L/D)  characteristics,  and  (3)  an  aero-configured  missile 
shape  having  "high"  (L/D)  characteristics.  Sketches  of  these  vehicles 
are  shown  in  Figures  C-l  through  C-3  and  the  aerodynamic  characteristics 
used  are  given  in  Appendix  B.  All  three  vehicles  were  used  in  developing 
data  set  1  and  the  medium  (L/D)  vehicle  was  used  to  develop  data  set  2. 

All  trajectories  were  initiated  at  300,000  feet  altitude  and  were 
terminated  at  a  speed  of  2,000  feet/second.  This  speed  termination 
allowed  the  terminal  altitude  to  vary  from  50,000  to  80,000  feet.  The 
data  shown  represent  aerodynamic  capability  only.  No  constraints,  such 
as  aerodynamic  heating  have  been  imposed. 

C-2 .  DATA  SET  1,  (L/D)  MAX  TRAJECTORIES 

Reentry  footprints  were  generated  using  (L/D)  max  and  the  single 
parameter  control  for  the  major  portion  of  the  footprint  boundary  and 
CL  max  for  minimum  down-range.  Initial  conditions  include  entry  velocities 
(inertial)  which  vary  from  20,000  to  24,800  feet/second  and  entry  path 
angles  from  -1.5  to  -25  degrees.  These  data  are  shown  in  Figures  C-4 
through  C-33.  Data  set  1  therefore  contains  a  complete  set  of  footprirts 
for  the  three  configurations  and  can  be  used  as  a  reference  for 
maneuvering  capability  available  for  a  range  of  configuration  types  and 
initial  conditions. 
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C-3 .  DATA  SET  2,  CONSTANT  ALTITUDE  TURNS 

This  data  set  is  for  trajectories  where  the  altitude  is  helc  constant 
during  the  turn.  All  trajectories  were  flown  at  zero  bank  angle  and  at 
L/D  max  until  the  initial  pull-out,  y  =  0.  For  straight  down-range  with 
no  turns,  the  angl e-of-attack  was  set  ’o  hold  level  flight  afte^  pull-out 
until  L/D  max  was  reached.  From  this  point  on  the  trajectories  were 
flown  at  L/D  max.  For  three-dimensional  trajectories,  turns  were 
initiated  at  the  pull-out  point.  The  angle-of-attack  was  set  tc  that 
required  for  L/D  max  and  level  flight  was  maintained  by  varying  the  bank 
angle.  The  footprints  were  generated  by  varying  the  amount  of  turn 
allowed.  Once  a  specified  heading  change  had  been  reached,  the  bank  was 
set  to  zero  and  the  altitude  allowed  to  vary  for  the  remainder  cf  the 
trajectory. 

Entry  conditions  include  variations  in  initial  velocity  (inertial) 
from  12,000  to  25,800  feet/second  and  entry  path  angles  from  -5  to  -25 
degrees.  Constant  velocity  (iso-velocity)  contours  are  shown  for  2000, 
4000,  6000,  8000,  and  10,000  feet/second. 

The  altitude  at  the  initial  pull-out  point  and  the  correspording 
range  are  shown  in  Figures  C-34  and  C-35,  respectively,  for  the  parametric 
variations  in  entry  velocity  and  path  angle  considered.  The  iso-velocity 
contours  are  shown  in  Figures  C-36  through  C-54.  The  dashed  lire  closing 
the  back  side  of  these  contours  represents  the  maximum  amount  of  turn 
possible  when  using  bank  angle  to  maintain  constant  altitude  du'ing  the 
turn  with  the  angle-of-attack  set  for  L/D  max.  To  obtain  a  faster  rate 
of  turn  than  shown  by  the  dashed  line  without  increasing  angle-of-attack 
a  steeper  bank  angle  would  be  required  allowing  altitude  to  decrease 
during  the  turn.  A  heading  can  be  obtained  for  any  of  the  possible 
trajectories  within  the  envelope  by  taking  a  tangency  point  to  this 
dashed  line  and  drawing  a  straight  line  through  the  iso-velocit.v  contours. 
An  example  is  shown  in  Figure  C-36  for  changes  in  heading  A a,  of  30,  45, 
and  90  degrees. 
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Data  set  2  (Figures  C-36  through  C-54)  then  contains  a  set  of 
footprints  for  the  medium  L/D  configuration  showing  the  maneuvering 
capability  available  if  the  turn  Is  performed  at  level  altitude  after 
pull-out.  The  pull-out  altitude  Is  the  same  as  used  in  data  set  1. 
Hence,  a  comparison  of  data  set  2  with  the  medium  L/D  configuration  of 
data  set  1  (Figures  C-14  to  C-23)  shows  the  large  loss  in  maneuvering 
due  to  the  constant  altitude  turn. 

For  example  comparing  an  entry  velocity  of  25800  feet/second  it  can 
be  seen  that  the  maximum  cross-range  of  3200  nmi  of  Figure  C-21  is  much 
larger  than  the  930  nmi  of  Figure  C-46. 
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Figure  C-l .  Low  (L/D)  Configuration 


C-2.  Medium  (L/D)  Vehicle 
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Figure  C-4.  Low  (L/D)  Footprint,  Entry  Angle 
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Figure  C-6.  Low  (L/D)  Footprint,  Entry  Angle 
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Figure  C-ll  .  Low  (L/D)  Footprint,  Entry  Angle 
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Figure  C-17.  Medium  (L/D)  Footprint,  Entry  Angle 
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Figure  C- 1 8 .  Medium  (l/D)  Footprint,  Entry  Angle 
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Figure  C-19.  Medium  (L/D)  Footprint,  Entry  Angle 
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Figure  C-2 5.  High  (L/0)  Footprint,  Entry  Angle 
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Figure  C-27.  High  (L/D)  Footprint,  Entry  Angle 
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DOWN-RANGE  ~  KN.MI. 

Figure  C-28.  High  (L/D)  Footprint,  Entry  Angle 
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Figure  C-30.  High  (l/UJ  Footprint,  Entry  Angle 
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Figure  C-33.  High  (L/D)  Footprint,  Entry  Angle 
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Figure  C-38.  Medium  (L/D)  Iso-Velocity  Contours,  (VT).=20,000  fps, 
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Figure  C-43.  Medium  (L/D)  Iso-Velocity  Contours,  (VT).=12,000  fps. 
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Figure  C-47.  Medium  ( L/D)  Iso-Velocity  Contours,  ( Vj )i =1 2,000  fps 
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Figure  C-52.  Medium  (I  /0)  Tso-Velocity  Contours,  (V.). =16, 000  fps. 


Figure  C-53.  Medium  (L/D)  Iso-Velocity  Contours,  (VT).=20,000  fps, 
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